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Abstract. - 

A layer of sand of thickness h flows down a rough surface if the inclination is larger than 
some threshold value Ostart which decreases with h. A tentative microscopic model for the 
dependence of Ostart with h is proposed for rigid frictional grains, based on the following 
hypothesis: (i) a horizontal layer of sand has some coordination z larger than a critical value 
Zc where mechanical stability is lost (ii) as the tilt angle is increased, the conflgurations visited 
present a growing proportion Za of sliding contacts. Instability with respect to flow occurs when 
z — Zs — Zc- This criterion leads to a prediction for 8start{h) in good agreement with empirical 
observations. 
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o 

A layer of dry sand on an inclined plane cannot remain stable if it is sufficiently tilted: 

there exists an angle Ostart above which grains must flow. Considering a grain on the top of 

►^ . the layer, this observation makes perfect intuitive sense j,: such a grain sits in a "hole" made 

^\ ' by the grains below, as the system is tilted the center of gravity of this grain shifts in the 

^^ . direction of the tilt, and the grain must eventually escape its hole, as illustrated in Fig.JTJa). 

Less intuitive is the observation 2 that Ostart decreases with the thickness h of the granular 

layer, and that this dependence can be observed for h as large as ten particle diameters, as 



C ' illustrated in Fig.([T]b). This implies that the stability of the layer is not governed by local 



rules: the grains at the top must be able to feel the presence of a fixed rough boundary 
significantly far below them. In general, one indeed expects that mechanical stability is a 
non-local property: imposing the stability of individual particles is a necessary condition, 
but is clearly not sufficient in general, as any collective motion of the particles must also be 
^^ , stable. The non-local microscopic criterion on the stability of an assembly of particles is due to 
^ ' Maxwell [5] . As we shall recall he showed that the key parameter is the coordination number 
z (the average number of contacts per particle) , which must satisfy some lower bound to allow 
stability. Here I shall extend this result to the situation where a fixed boundary (the rough 
surface of the inclined plane) and a free boundary (the upper layer of grains) are present. This 
result, together with the hypothesis supported by numerical analysis 4 that more and more 
contacts are sliding as the inclination increased, leads to a reasonable estimate of 9start{h)- 

I start by recalling Maxwell's criterion for the mechanical stability of N frictionless particles 
in d dimensions, forming Nc contacts, whose coordination is z = 2Nc/N. If kij is the stiffness 
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Fig. 1 - (a) A disc sitting on a regular array of similar discs is unstable if the tilt angle is larger than 
dstart = 30 . This simple picture cannot explain the dependence of Bstart with the thickness of the 
granular layer, (b) Layer thickness in units of particle diameter h/do v.s. 9 start in an incline plane 
configuration where grains flow down a rough surface. Crosses correspond to the empirical results of 
Forterre and Pouliquen [5] . The curve is the fit of the prediction Eq. (|12p . 



interaction, and if the strain of the contacts is small [6j , the energy expansion can be written 
as: 



k 



SE = J2 -^[{SR. - SRj) ■ nyf + oiSR^) 



(1) 



iv) 



where n^ is the unit vector going from i to j, and SRi is the displacement of particles i. A 
system is mechanically stable if it cannot be deformed without an energy cost, i.e. there exists 
no soft mode or particle displacement for which 6E ~ 0. Using Eq.(IT]) this implies that the 
system of equation 

{SRi-SRj)-n,j ^Oyij. (2) 

must not present any non-trivial solution. Eq.([2]) is a linear system and non-trivial solutions 
must exist if it has more degrees of freedom than equations. There are Nd degrees of freedom 
and Nc equations, therefore stability requires Nc > Nd or equivalently z > 2d, as derived 
by Maxwell. Thus, any solid made of frictionless particles must satisfy z > 2d. It turns out 
that in the limit of large stiffness and fixed compression, the system must become isostatic, 
i.e 2; = 2d [H [II 1^. This can be understood as follows: for infinite stiffness particles are just 
touching each other, therefore each pair of particles in contact must satisfy: 



\R.,-Rj\\ = doyij. 



(3) 



where do is the particle diameter. Again Eq.Q has Nd degrees of freedom and Nc constraint, 
it can therefore be satisfied only if Nd > Nc. Together with the stability condition Nc > Nd, 
this leads to Nd = Nc oi z = 2d, as observed numerically [TUj. Thus frictionless rigid grains 
have a well-defined coordination. 

For frictional particles the situation is different. As far as stability is concerned, there 
are now d{d + l)/2 degrees of freedom per particles (d translations and d{d — l)/2 rotations), 
and each contact now brings d constraints ( 1 due to the radial interaction and d — 1 due 
to friction, corresponding to the directions transverse to the contact). The Maxwell criterion 
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reads z > d+ 1 = Zc- In the stiff limit, the argument associated with Eq.Q is stiU vahd, and 
therefore we get for the coordination: 

d+1 <z <2d (4) 

In frictional rigid grains the coordination is therefore not fixed, but depends on the system 
preparation, and on tlie friction coefficient /i. Numerically [HI [TT] it is found that as /i ^ 
leads z ^ 2d (i.e when the friction coefficient vanishes one recovers the frictionless limit), 
whereas /x — > oo leads to z — > d + 1. For example for poly-disperse discs with /i = 0.5 one 
finds numerically z « 3.2 [TT] . 

In a frictional assembly of grains some contacts can slip: the ratio of transverse to lon- 
gitudinal force in such contacts is fi. These contacts display no more transverse stiffness, 
they do not bring transverse restoring force when an external load is applied. If Ns is the 
number of slipping contacts, the system displays Ns{d — 1) less constraints than if no contacts 
were slipping. We shall assume that slipping contacts leave the number of degrees of free- 
dom relevant for the mechanical stability unchanged (in the discussion section we shall argue 
that this assumption is valid for intermediate and large friction coefficient, but not for nearly 
frictionless particles). Introducing: 

Zs = 2id-1)N,/N, (5) 

the Maxwell counting arguments now leads to: 

Z> Zc + Zs (6) 

Under increasing shear, or when a granular layer is tilted, one expects for the configurations 
visited that the ratio of the typical transverse contact force to normal force will grow. This 
leads to an increase in the number of contacts that slip. In the limit of rigid grains the 
dependence of Zg on the applied stress can only be a function of the ratio a/p = tan{9), where 
p is the pressure and a the shear stress: 

z, = /^(ton(0)) (7) 

Numerical computation of f^itaniO)) in two dimensions for jjl = 0.5 are available[4]. Eq.® 
and Eq. ([7]) leads to the following condition on 6 for the stability of an infinitely thick layer of 
grains {h = oo): 

f^{tan{0)) <z-z, (8) 

The coordination z may a priori depends on 9. We shall assume that this is not the case, 
as supported by numerics [4 . Relaxing this hypothesis would lead to no qualitative changes 
in the following. We also assume that a layer of sand starts flowing when the bound® is 
saturated: 

ff_^{tan{9 start{h = oo))) !V Z - Zc (9) 

This assumes that the destabilization induced by slipping contacts occurs at an angle where 
other destabilizing effects (such as the local instability described in Fig.([TJa)) did not yet set 
in. As far as orders of magnitude are concerned our assumption is supported by the numerical 
studies [4l[TT] observing respectively f^i{tan{9 start)) ~ 0.16 and z — Zc ^ 0.2 for polydisperse 
discs of friction ^ = 0.5. In other words, Eq.® leads to a reasonable estimate of 9start in 
those numerics. 

We know extend this result to the case of finite thickness, and first show that the presence 
of a fixed boundary (such as the rough plane on which grain sits) increases the coordination 
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by some amount ado/h, where h is the thickness of system and a is a constant of order one. 
For concreteness consider a system of hnear size h of N frictionless particles with periodic 
boundary conditions (the argument also holds with friction for any geometry), such as the 
one shown in FigOa. If m disjoined particles are pinned, the system formed by the other 
particles loses md degrees of freedom, but still present the same number of contacts. The 
ratio of constraints per particle has therefore increased: this is equivalent to increasing the 
coordination by an amount Az = 2[Nc/{N - m) - Nc/N] w 2{Nc/N){m/N) = zm/N, where 
we have assumed m/N << 1. To mimic the presence of a rough surface we can pin all the 
particles crossing a surface of the system, such as the red line represented in Figj^lb. In this 
case a fraction of order m ^ dg/h of the particles are pinned, and the coordination increases 
by some amount Az = ad^/h with respect to the coordination z of a bulk solid. 

To study the effect of a free boundary we shall instead consider the situation where m 
particles are removed of the system. The remaining system looses again md degrees of freedom, 
but it also looses mz contacts. As far as degrees of freedom and constraints are considered, 
this is equivalent to decrease the coordination by some amount Az = 2[{Nc — mz)/ {N — m) — 
Nc/N] « —am/N. To mimic the presence of a free boundary, we can again draw a surface and 
remove all the particles crossing it, this procedure creates a free surface, as shown in Figlllc. 
From this simple argument we may (too) naively expect that the effect of a free boundary 
is exactly the opposite a fixed boundary: it decreases the coordination by Az = —ada/h. If 
so, the finite size effect of a layer of grains with a free and a fixed boundary would simply 
cancel out. This is nevertheless not the case, because as it appears in Fig[2lc, removing a 
layer of particles does not lead to a realistic description of a free interface: in particular a 
finite fraction of the particles are left with less than d contacts, an unrealistic situation. A 
better description of a free interface is obtained by letting such particles flow down, until 
they make d contacts, as exemplified in Fig[21d. Since this occurs for a finite fraction of the 
particles at the interface, the coordination of the system increases by some amount of order 
do/h. Consequently the presence of a free boundary decreases the effective coordination by 
an amount Az = a'da/h, with a' < a. 

Accounting for these effects, Eq.® can be generalized to describe a layer of grains with 
finite thickness: 

ff^{tan{0 start)) = z- Zc + {a- a')da/h (10) 

It is convenient to expand f^{tan{9)) in the vicinity of dstart{h — oo) where the value of /^ is 
known from Eq.Q: 

f^{tan{9)) « (z - Zc) + a^{tan{9) - tan{9 starti'x,))) (11) 

where a^ is the derivative of /^^ at 9start(po). a^ is expected to be a decreasing function of /i, 
as a larger friction coefficient will decrease the rate at which the number of slipping contacts 
increases. Using Eq. (fTT|) . Eq. ([TO| can be rewritten as: 



do a^{tan{9start{h)) - tan{9 start{oo))) 



(12) 



which is our final result. Once the value of 9start{oo) is fixed (and as discussed above, Eq.Q 
leads to a reasonable estimate of this quantity in the numerical systems where it can be tested) , 
Eq. (|12p has one remaining fitting parameter (a — a')/a^, which we expect to be an increasing 
function of the friction coefficient. Fig[T]b shows that the form proposed in Eq. p^ fits well 
the empirical observations. 

I conclude with several remarks: 
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Fig. 2 - (a) Assembly of frictionless particles with periodic boundary conditions (b) Model of a rough 
substrate: all the particles intersecting an arbitrary horizontal line, represented in red, are pinned. 
This increases the stability of the system, in a manner equivalent to an increase of coordination, as 
computed in the text, (c) Naive model of a free interface, generated by removing all the red particles 
crossing the arbitrary line. This model is not realistic because some of the particle are left with 
less than d contacts, and are therefore unstable. The presence of such an interface decreases the 
stability of the system, and is equivalent to a decrease of coordination. The effect is exactly opposite 
to the presence of a fixed boundary, (d) A more realistic free interface is generated by letting the 
unstable particles at the surface flow and make d contacts. This increases the coordination. As a 
consequence the destabilizing effect of a free boundary is milder than what expected from the naive 
model represented in (c), and the overall effect of the joined presence of a free and fixed boundary is 
stabilizing. 



(i) The model proposed here may apply to a rough plane, such as the one obtained by 
gluing particles on a surface. If the interaction grain-plane is instead described by a friction 
coefficient /i^ significantly smaller than /x, one expects that the first contacts to slide will be 
at the interface grain-surface, and that the granular layer will slip as a bloc for some angle 
determined by /z^, and Ostart should not depend significantly on h. At the other extreme 
where the plane is extremely rough, the grains in contact with the plane may be significantly 
immobilized. If so the present result can still be applied replacing hhy h — do. 

(ii) Our analysis neglects spatial heterogeneities in the coordination and in the location 
where contacts slip, and is in some sense a mean-field model. If large fluctuations in these 
quantities occur, some regions may yield for angles smaller than our estimate. Small avalanches 
are indeed observed as the inclination increases, and become more frequent close to Ostart |13j . 
Our analysis will hold qualitatively if such avalanches cannot trigger a global flow except near 
the critical inclination for which the system as a whole becomes unstable. 

(iii) In a shear cell where the two boundaries are fixed, the same line of thought leads to 
an equation similar to Eq. (|12p . where a — a' is replaced by 2a. One thus expects the finite 
size effects to be stronger in this geometry. 

(iv) When deriving Eq.®, we have assumed that the contacts that slip leave the number of 
degrees of freedom relevant for mechanical stability unchanged. This is true except when the 
fraction i/ = ff_i/z of such contacts is large. In that case, all the contacts of a particle may slip, 
with a probability of order u^ . Rotating such a particle leaves the contact forces unchanged. 
Thus as far as linear response is concerned these particles act as frictionless particles, and their 
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degrees of freedom of rotation is irrelevant for mechanical stability. This implies that Eq. ^ 
over-estimates the destabilizing effects of the slipping contacts. Nevertheless for intermediate 
and large friction coefficient, the corresponding corrections are negligible. For example in 
numerics j4j with fi = 0.5 it was observed that i/ « 0.05 and z ~ 3, and we expect the number 
of particles acting as frictionless ones to be of order 0.05^ = 1.25 * 10^^ << 1. Nevertheless, 
this effect must be dominant if frictionless particles, such as emulsions, are considered. Then 
one expects the system to remain isostatic (z = 2d) independently of 0, leading to /^ — 0. 
For emulsions the nature of the instability at dstart rnust therefore differ from the present 
description. 

(iv) Once the flow is initiated, it stops when the slope of the granular surface becomes 
smaller than some angle ^^top < Q start- Empirically Ogtop depends on /i in a way qualitatively 
similar to 9 start p] , and it seems unlikely that these two similar dependences stem from very 
different causes. This suggests that the properties of dense granular flows are affected by the 
presence of nearly stable configurations, an idea with long history in the context of the glass 
transition [M]. 

It is a pleasure to thank B. Andreotti, A. Kabla and O. Pouliquen for discussions and 
comments on the manuscript. 
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